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7-H , Abstract 

We consider nonholonomic geodesic flows of left-invariant metrics and 
<*i . left-invariant nonintegrable distributions on compact connected Lie groups. 

Oh' The equations of geodesic flows are reduced to the Euler-Poincare-Suslov 

equations on the corresponding Lie algebras. The Poisson and symplectic 
structures give raise to various algebraic constructions of the integrable 
Hamiltonian systems. On the other hand, nonholonomic systems are not 
Hamiltonian and the integration methods for nonholonomic systems are 
much less developed. In this paper, using chains of subalgebras, we give 
constructions that lead to a large set of first integrals and to integrable 
K^ , cases of the Euler-Poincare-Suslov equations. Further, we give exam- 

$_^ ■ pies of nonholonomic geodesic flows that can be seen as a restrictions of 

integrable sub-Riemannian geodesic flows. 
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Introduction 

In this paper we are interested in the geometry and integrability of nonholo- 
nomic geodesic flows of left-invariant metrics and left-invariant nonintegrable 
distributions on compact Lie groups. The equations of geodesic flows are re- 
duced to the Euler-Poincare-Suslov equations on the corresponding Lie alge- 
bras. These systems are natural generalizations to Lie algebras of the Suslov 
nonholonomic rigid body problem and were introduced independently by Ko- 
zlov IJJl and Koiller [ ftl|. The known integrable cases are given by Fedorov and 
Kozlov H and author |ifj. 

In the recent years appears many papers concerning geometrical formula- 
tion of the nonholonomic mechanics. For instance, see pi pi O] and references 
therein. However, since nonholonomic systems do not admit a Poisson struc- 
ture, the integration theory of the constrained mechanical systems is much less 
developed then for the unconstrained. We mention Chaplygin's results M. By 
the use of an invariant measure, he had given some of the most interesting ex- 
amples of the solvable nonholonomic systems and noticed that the phase space 
could be foliated on invariant tori, placing these systems together with inte- 
grable Hamiltonian systems. The methods of integration of systems with an 
invariant measure, as well as illustrative integrable examples can be found in 

ill ill- 

Now, we shall briefly describe the results and outline of this paper. 

In sections 1 and 2 we shall give basic definitions and notation. 

One of the well known ways for studying the integrability of Riemannian 
geodesic flows is by using certain nitrations of Lie algebras (sec ^, [l4|, [L6| |l H) . 
By taking appropriate chains of subalgebras, we can get integrable cases of the 
Euler-Poincare-Suslov equations. 

In section 3 we shall consider the case when the left-invariant nonintegrable 
distribution is an invariant subspace of left-invariant metric. Then reduced 
system has an invariant measure. We shall construct integrable examples by 
using chains of the form K C H C G, where (G, H) is a symmetric pair. The 
obtained systems are generalizations of the Fedorov-Kozlov integrable case. 

In section 4 we deal with an arbitrary chain Go C G\ C . . . C G„ = G 
This gives us opportunity to construct exactly solvable examples without an 
invariant measure as well. 

In some cases the interesting phenomena arises: the nonholonomic geodesic 
flow could be seen as a restriction of Hamiltonian flow on the whole, uncon- 
strained phase space. We shall present here two families of integrable sub- 
Riemannian geodesic flows which restrict to the our non-Hamiltonian problem 
(section 6). 

1 Nonholonomic geodesic flows 

Let (Q, (•, •)) be n-dimensional Riemannian manifold with Levi-Civita connec- 
tion V. Let T> be the nonintegrable p-dimensional distribution distribution of 



the tangent bundle. The distribution can be defined by n — p independent one 
forms o,i in the following way: 

V q = {£eT q Q, ai(O=0, i = l,...,p}. 

The smooth path j(t),t 6 A is called admissible (or allowed by constraints) 
if velocity j(t) belongs to £> 7 (t) for all i e A. There are two approaches to 
define the geodesic lines among admissible paths: by induced connection as a 
" straightest" lines and by variational principle as a "shortest" lines. We shall 
deal with the first approach which arises from mechanics. The admissible path 
7(i) is called a nonholonomic geodesic line if it satisfied d'Alambert-Lagrange 
equations: 

*r(V.y (t) 7(t)) - 0, (1) 

where it : T q Q — > V q , q G Q is the orthogonal projection. 

Let {•, •} be canonical Poisson brackets on T*Q. By Xf we shall denote the 
Hamiltonian vector field of the function / (dg(Xf) = {g, /}, for all g : T*Q — > 
R). Let M. be the constraint submanifold in the phase space T*Q: 

M = {{p, q) G T*Q, p G g q (V q ) C T^Q}, 

where we consider scalar product (•, -) q as the mapping g q : T q Q — » T*Q. Taking 
for the Hamiltonian h(p,q) = ^p(g,7 1 p), p G T*Q, we can write (pj) in the 

following form: 

p 

x = X h (x) + y^ Aivertai| .( x ) , (2) 

where verta^^) G T X (T*Q) is "vertical" lift of a,| CT ( x ) G T*,-,Q, a : T*Q -» 
Q is the natural projection and Lagrange multipliers Ai are chosen such that 
the phase trajectory x(t) belongs to A4. In canonical coordinates x = (p,q) 
equations (g) are: 

9h{p, q) y~y . dh(p, q) 

It is important to note that the equations are not Hamiltonian (they are 
Hamiltonian with respect to the almost-Poisson brackets on M |L3L g]) and 
that the Hamiltonian function is conserved. 

Sufficient condition for the integration of the non-Hamiltonian system 

x = f(x), xel™ (4) 

is existence of to — 1 independent integrals, or one integral less in the case of 
the existence of an invariant measure. 

Suppose that the system of equations (0) has an invariant measure and 
771 — 2 first integrals F\ 1 . . . , F m . If Fi, . . . , F m are independent on the invariant 
set M c = {x G M m , Fi(x) = a, i = l,...,m - 2} then the solution of (g) 



lying on M c can be found by quadratures (the Jacobi theorem). Moreover, if 
L c is a compact connected component of M c and f(x) ^ on L c then L c is 
diffeomorphic to a two-torus; one can find angular coordinates (pi, </?2 on L c in 
which equations (0) take the form similar as in the Liouville theorem: 



Wl . U)2 

Pi = ^7 \> ^2 



$(</?!, ^2)' ~ $((/?i,yJ 2 )' 

where wi , uji are constant and $ is a smooth positive 2-7r-periodic function in 
tpi,tp2 (see 0). 

Therefore it is natural to call a non-Hamiltonian system integrable if it can 
be integrated by the above procedure; or more generally (as it was point out in 
|{22|), if the trajectories of the system belong to invariant tori with dynamic of 
the form 

The flow of (ra) is unevenly winding and admit the invariant measure fi(D) — 
J D <&d<pi A ... A difik- Note that it is shown in (fl that for almost all frequencies 
Wi, . . . , u>k, by smooth change of variables 

T fc {^,...,^}^T fc {0 1; ...,0 fc }, 

equations (ph can be reduced to the form 6i = flj = Wj/II, i = 1, . . . , k, where 
II denotes the total measure of T fc . 



2 Euler— Poincare— Suslov equations 

Now, let Q be a compact connected Lie group <S with Lie algebra G = T e <8. In 
what follows we shall identify G and G* by Ad® invariant scalar product (•,•}; 
T© and T*© by bi-invariant metric on ©. 

We shall consider left-invariant distributions. Let 

D = {lueG, {lu, a 1 ) = 0, i = 1, . . . , p} C G 

be the restriction of the left-invariant distribution V to the Lie algebra G, for 
some constant, linearly independent vectors a % in G. The distribution is nonin- 
tegrable if and only if D is not a subalgebra. From the invariance, we have that 
■D g = g ■ D. 

Let I : G — > G be a symmetric, positive definite (with respect to (•,•)) 
operator that induces left-invariant metric: 

(m,V2)g = (I(u>x),u>2), r\\,r\2 G T g &, u>i,u>2 G G, r]i= g-u>i, r) 2 =g ■ u> 2 - 

Let M be the restriction of the constraint submanifold M. to G. Then A/ = 
/(D) and .Mg = g ■ M. The Hamiltonian of the geodesic flow is function 
h : TG5 — -> R obtained from reduced Hamiltonian function H{x) = ^(A(x), x) by 
left-translations (here A = I -1 ). 



In a such notation, the equations (0) are reduced to: 

p p 

x = [x, VH{x)] + Y^ >*<? = [«, A (x)] + Yl X ^' ( 6 ) 

j=l i=l 

where Lagrange multipliers are chosen such that x belongs to M = 1(D), i.e. 
such that lo — A(x) belongs to D: 

{A(x),a l )=0, i=l,...,p. 

According to Fedorov and Kozlov || we shall call these equations the Euler- 
Poincare-Suslov equations, as a generalization of the Suslov nonholonomic rigid 
body problem. They have a quite different nature of the corresponding Euler- 
Poincare equations x = [x,A(x)Y. For instance, in general, they do not have a 
smooth invariant measure (see [H). The equations (ph could be seen also as a 
reduced equations from a point of view of a reduction of nonholonomic systems 
with symmetries given in g, [L|] . 

Nonholonomic geodesic lines g(t) are solution of the kinematic equation 
g~ x {t) ■ g{t) — uj(t) — A(x(t)), where x(t) are solutions of (pi). In other words, 
the following diagram commutes: 

V t 

M -> M 

A | | A (7) 

M -» M 
p t 

Here, V 1 and P t are phase flows of the nonholonomic geodesic flow and Eulcr- 
Poincare-Suslov equations; A mapsp = g-x G T g <3 to x G G. If Euler-Poincare- 
Suslov equations are integrable we shall say that nonholonomic geodesic flow is 
integrable in the sense of the factorization (0). 

The reduced Hamiltonian function H(x) = \{x,A(x)) is the first integral 
of system (o). This follows from the conservation of energy in the natural 
mechanical nonholonomic systems with linear constraints. One may also prove 
this fact by a direct computation of the time derivative of H(x) along the 
vector field defined by (ph. Namely, integral F of Euler-Poincare equations 
x = [x, A(x)} is the integral of Euler-Poincare-Suslov equations (m) if and only 
if: __ 

J2M^F{x),a t )\ xeM = 0. (8) 



3 Symmetric pairs 

Let L be the subspace of G spanned by a 1 , i = 1, . . . , p. In this section we shall 
consider the case when A preserve the orthogonal decomposition G = L + D, 



i.e., A = Al + Ad, where Al : L — » L, Ad : D — > D arc positive definite 
operators. Then M = 1(D) = D and we can write (@) in the following way: 

£=[t,A D (£)] D , £eD (9) 

(by xk we denote the orthogonal projection of x to the linear space K). 

The equations (g) preserve the standard measure on D. Also the constrained 
reduced Hamiltonian function Hr> = -^ (£, .Ad(£)) and invariant -F(£) = (£, £) are 
always first integrals of the system. Note that by (g), in general, the invariant 
F(x) = (x,x) is not the integral of (|6[). 

Example 3.1 We have, at least, following integrable cases: 

1. If Ad = s ■ I do, s£l then the solution of (0) are £ = const. In this case 
the constraints have no influence to the motion (the Lagrange multipliers in 
(ra) are equal to zero). The nonholonomic geodesic lines are simply given with 
g(t) = goexp(nt), r) € D, go € 0. At the same time these lines are geodesic 
lines of the left-invariant metric induced by A = A^ + s ■ I dp . From now on we 
shall suppose that Arj ^ s ■ Idrj . 

2. dimD = 2. Then the solution of (||) are £ = const. 

3. dimD = 3,4. Then the solutions of (||) belong to the intersections of the 
spheres and ellipsoids: 

M C1 , C2 = {£ed, h d (0 = ±(s,a d (0) = a, F(0 = (tO = c 2 }. 

For dimD = 3 general trajectories are periodic. If dimD = 4 then (||) could 
be integrated by Jacobi theorem as a system with an invariant measure (note 
that general connected components of invariant submanifolds M CliC2 are 2- 
dimensional spheres). 

4. Recall that (G, H) is called a symmetric pair if the following condition is 
satisfied: 

[H,H}CH, [H,V]CV, [V,V]CH, 

where V is the orthogonal complement of H . If there is a subalgebra H c G 
such that (G,H) is a symmetric pair and that H C L C G, then [D.D] C H 
and [D,D]d = 0. Therefore the solution of (O) are ^ = const. 

In the cases 2 and 4 above, nonholonomic geodesic lines are g(t) = go exp(?7i), 
r\ £ D, go £ (3, but these lines not need to be geodesic lines of the left-invariant 
metric induced by A. 

Motivated with the last example, let us consider the chain of subalgebras: 

KcH CG, 

such that (G, H) is a symmetric pair and that H is not a subspace of L. Let 
G = H + V be the orthogonal decomposition. We can consider the adjoint rep- 
resentation of K on the linear space V: ijeif^ [r/, •] e End(V). Decompose 1/ 



on irreducible subspaces V = Vq + V\ + . . . + V m . Here Vq denotes the subspace 
with trivial representation. 
Let D be of the form: 

D = U + W + W 1 + ... + W rn , U = DnK, W k = DnV k . 

Since H <£. L we have that dimC/ > 1. Let dimWfc > 1, k = l,...,n, 
dim Wk = 0, k > n. 

Theorem 3.1 Suppose that operator Ajy preserve the decomposition D = U + 
Wq + W\ + . . . + W n and that Ar>\u = s ■ Id\j, S G R. Then equations ([£), 
besides functions -ffo(£) = \ ( A £> (£),£) and F(£) = (£,£)> have a set of the first 
integrals, the projection of £ to Wq : Fq (£) = £w an d functions: 

Fk(0 = (Bw k (tw k ),Zw k ), k=l,...,n 

where B Wk — A D \ Wk - s ■ Id Wk . 

Proof. Let A w — A D \ W , A Wk = A D \ Wk , W = W + . . . + W n . The equations 
(^|) have the form: 

j t {£u + £w) = [tv + £w,8tu + A w (tw)} D . (10) 

From [17, V] C V, [W, W] C H and @ we get: 

£,u = [£wiAw{£w)]u, 

£w = [£u,Aw(t;w) — st;w]w- 

Since Aw preserve the decomposition W = W^ + W\ + . . . + W ni the second 
equation is separated on n + 1 equations: 

^w = 0, £w k = [Zu,Bw h (€w k )]w k , k=l,...,n. (11) 

It is clear that Fk are integrals of (O). Note that the invariant F = (£,£) is 
dependent of functions Ho and F^, k = 0, . . . , n. 

Corollary 3.1 If operators B\v k are positive definite and c„+i satisfied inequal- 
ity: 

n 

Cn+i > |co| 2 + y"T-, b k = min {B Wk (€w h ),€w k ), 
k=1 b k \tw k \=l 

then invariant subspaces: 

M c = {£eD, $ Wo = c , Fi(0 = Cl) . . . , F B (0 = c„, F(0 = c„ +1 } 
are diffeomorphic to the product of spheres: 

odim Wi — 1 v v CfdimWn— 1 v odinit/ — 1 

In particular, if diuiWk <2,k=l,...,n, dim C7 = 1, i/ien M c is diffeomorphic 
to the disjoint union of g < n dimensional tori with quasi-periodic dynamic (fq,). 



Proof. The first part of the corollary follows from the relations: 



|2 „ |„ |2 V^ i e |2 it |2 s- ° k 

- 2^ \€w k \ 



ICul 2 =Cn+l - |c | 



fc=l 



£w fc | < -r- , (Bw k (€w K ) ) £,w k ) =Ck- 

Ok 



Let dim Wk = 2, k = 1, . . . , g, dim Wk = 1, fc = g + 1, . . . , n. Let </>& mod 2tt 
be the angular variables of ellipses Fk — Ck, k = 1, . . . ,g. Let £[/ = $ ■ 77, $ e M, 
7/ G C/, |?/| = 1. Then $ can be expressed in terms of the <pk up to the sign: 



*-*(0,,...,*,)-± 



\ 



c„+i - |co|- -2^l ?ir - ■ l " A 
fc=i 



5> 



fe=9+l 



The sign of $ is determined by the choice of the connected component T 9 
of M c . 

By substitution of £iy fc = £w fe {<f>k), k = 1, . . . ,g, & = $(<£i, ... ,<j) g ) ■ rj to 
(y_l|), we obtain that angular variables satisfied the following type equations: 



ii = *( 



;i, . . . , (p g 



0fl)/l(0l).....0ff = *( 



<t>g)fg{<i>g)- 



>l,.-.,9g)Jg\9g 



(12) 



The system ( |12j ) can be integrated in term of the new time r defined by: rfr = 

If all /fe are different from zero, then we can introduce angular variables ipk 
by averaging: 



Vfe = <i?k{4>k) = u k 



ds 



/o A-(s) 
in which ( |12|) takes the required form 



w fe = 2ir 



^ ds - ' 



fk(s) 



, k = l,...,g, 



<Pi 



$- 



.Vs 



$- 



(13) 



Remark 3.1 The frequencies Wj depend only of the metric A. If the trajectories 
are periodic on one torus, they are periodic on the rest of the tori as well. 



Remark 3.2 We have conditions dimWi = dimVi = 2, dimU = dimK = 1 
taking: 



G = 



u wt W L ). r = .s„(2>. /. = .s,x/;-2i. D = r + n: 



So, we can see the corollary as a generalization of the Fedorov-Kozlov integrable 
case [pi. 



Remark 3.3 Besides of preserving the standard measure in D, we can easily 
rewrite system (ph in the " Hamiltonian form" : 

F = {F,H D } D , F:D->R, (14) 

where Hd = \{£,, Ad(£)) an d {■, -}d are almost-Poisson brackets defined by: 

{F ll F 2 } D (0 = (^lVF 2 (t),VF 1 (0}), F U F 2 :D^R. (15) 

These brackets are bi-linear, skew-symmetric and satisfy the Leibniz rule. In 
the general case, they do not satisfy the Jacobi identity. For D = G these are 
the usual Lie-Poisson brackets on the Lie algebra G. 

As for the Poisson brackets, we can define central (or Casimir) functions of 
the brackets (|15|). These are the functions F that commute with all functions 
F : D — > K. Obviously, they are integrals of the system (^). The example of 
the central function is F(£) = (£,£). 

More about the almost-Poisson setting for the nonholonomic systems can be 
found in || §. 

4 Chains of subalgebras 

Suppose we are given a chain of connected compact subgroups: 

©o C 0i C . . . C <8„ = 

and the corresponding chain of subalgebras: 

Gp C Gi C . . . C G n = G. 

Let Gi — Gi-i + Vi be the orthogonal decompositions. Then 

G l = G + V 1 + ... + V l . 

Following ||, Q , consider A of the form: 

A = A + si • Id Vl + ... + s„ ■ Id Vn , s t > 0, i = l,...,n, (16) 

where Aq is a symmetric positive operator defined in the subalgebra Go- 
Suppose that D has orthogonal decomposition: 

D = D + D 1 + ...+D n , (17) 

D k = {^k e 14, (a\,ojk) = 0, i = 1, .. .,Pk}- 

Then 13^, fc > are invariant subspaces of A. By x k denote the orthogonal 
projection of x to Dk, fc > 0; and by xq denote the orthogonal projection to Gq. 



Theorem 4.1 The Euler-Poincare-Suslov equations ^), with D and operator 
A(x) of the form j\lq ) and fcP\), are equivalent to the Euler-Poincare-Suslov 
equations on the Lie subalgebra Gq: 



Pa 



%o = [x ,Ao(x )] + > mob, (18) 



i=i 



(A (x ),a l ) = 0, i = l,...,p , 

together with a chain of linear differential equations on the subspaces D k : 

x k = [xk, A Q (x ) - s k x Q + (si - s k )xi + . . . + (sfe-i - Sk)xk-i]k- (19) 

Proof. Let xy k be the orthogonal projection of a; to 14. The simple compu- 
tations show: 

[x, A(x)] — [x + x Vl + . . . + xv„, Aq{xo) + SiX Vl + . . . + s n xv„] 

= [x Q ,A (x )} + [x Vl ,A (x ) - siXq] + .. .+ (20) 

+[xv n ,A (x ) - s n x + (si - s n )x Vl + . . . + (s„_i - s„)xy„_ 1 ]. 



Taking into account ( fL7| ) , 
orthogonal projection of 
follows. 



and relations [G k -i, V k ] C 14, k = 1, . . . ,n, after 
to the linear space Go + D\ + . . . + D n the theorem 



Remark 4.1 Let D k = 14 and let T k be the algebra of polynomials that are 
constant on the orbits of the adjoint action of the group (5fe_i on 14- It can be 
easily seen that they are integrals of (119) . The number of functionally indepen- 
dent polynomials in J- k is equal to: 

dim 14 - dimGfc-i + min dim{y <E G k -i, [y,Xk-i] = 0}. 

Xk&V k 



If Eulcr Poincare-Suslov equations (|18|) on Go are solvable then the inte- 
gration of original equations rfq) reduced to successive integration of the chain 
of linear dynamical systems ( J19| ) for k > 0. 

The most simplest case is when the solutions of ( |lq ) are xo — const. Then 
the vector x\ satisfied a linear equation with constant coefficient and it is a 
elementary functions of the time t. This is happened if Aq — Ida or if Go is a 
commutative subalgebra (see also example 3.1) In particular if dimZ?o = then 
we have that xq = 0. In that case x,\ = and X2 is a elementary functions of 
the time t. This case will be treated again in the last section (example 6.1). 

Example 4.1 In addition, we can obtain exactly solvable cases without an 
invariant measure. For instance, take the chain: 

so(3) C so(4) . . . C so(n). 

Let Dq — {luq g so(3), {ujq, oq) = 0}. If ao is not an eigenvector of the Aq then 
equations: 

x = [xq, A (x Q )} + Xa Q , (A(x ),a ) = 0, (21) 



10 



have no invariant measure JL2[ . After identification (so(3), [•,•]) = (R 3 , x), the 
system pl| ) describe the rotation of a rigid body fixed at a point and subject 
to the constraint: the angular velocity ujq is orthogonal to the fixed vector in 
body coordinates do. This nonholonomic problem was solved by Suslov Jl8[ . 

The phase space Mq = Iq(Dq) has the following property: there is a asymp- 
totic line I such that lim^-too Pq(SIo) C I, where Pg is a phase flow on Mq and 
flo is any subset of M . Then the phase flow P t of f[§|), ( |i9| ) has asymptotic 
hyper-plane and also has no invariant measure. 

However, for dimD/j < 2 we can easily integrated corresponding systems 
(P^l). Suppose that we know yk-\{t) = A (:eo(£)) - s k x (t) + (si - s k )xi(t) + 
. . . + (sfc_i — Sfe)xfc_i(i) as a function of time. If dimPfe = 1 then x k — const. 
For dimPfe = 2, let 0^ be the angular variable of the circle I k = (xi^Xk) = c&. 
Then it can easily be checked that equations x k = [xk,yk-i(t))k get the form 
4>k = fk(t), where fk{t) is some known function of time. Note that connected 
components of invariant submanifolds H = h,I\ = Ci,...,I n — 3 = c n ~3 are 
diffeomorphic to tori, but the dynamic on tori is quite different from (J13J) . 

5 Reconstruction 

The integrability of Euler-Poincare equations x = [x, A(x)] implies non-commu- 
tative integrability of unconstrained geodesic flow — the phase space T<3 is 
foliated on d < dim 05 dimensional invariant isotropic tori with quasi-periodic 
dynamics (see flq]). In the nonholonomic case there is no Poisson structure. 
In order to precisely describe dynamic on the whole phase space M. we have 
to solve kinematic equation g~ l {t) ■ g(t) = ui(t) = A(x(t)). This problem, for 
the Fedorov-Kozlov integrable case, is studied by Zenkov and Bloch |24|, |25f. 
Similar analyses can be applied to the considered integrable cases: 

(i) The periodic solutions on M correspond to quasi-periodic motions on the 
phase space M. 

(ii) Suppose that in the corollary 3.1 we have Ck = eCfc k = 1, . . . , g, where 
e is a small parameter, and that the following Diophantine conditions hold: 

\i + ifam + ... + k g u; g )\ > cflfcii + . . . + iyr 7 , 

1 = 0,1,2, (k 1 ,...,k g )eZ9-{0} 1 

for some constants c > and 7 > g — 1. Then the reconstruction of quasi- 
periodic motion ( |l3|) to T> = M can be approximated by quasi-periodic dynam- 
ics on the time interval of length ~ exp(l/e). 

6 Hamiltonian cases 

In some cases, the nonholonomic geodesic flow (0) on M can be gotten as a 
restriction of Hamiltonian flow from T*Q to submanifold M.. We shall say that 
nonholonomic geodesic flow (0) has a Hamiltonian restriction property if there 
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is a function h* : T*Q — > M such that the restriction of Hamiltonian equation 
x = Xh* to M coincide with (||). 

The basic examples are when Lagrange multipliers in (0) vanish. Then we 
can take h*(p,q) = h(p,q) — \p(gq X p)- However, we shall see that there is 
another natural choice of the Hamiltonian h* . Note that the cases 2 and 4 of 
example 3.1 have Hamiltonian restriction property with Lagrange multipliers 
(in general) different from zero. 

From now on, we shall use the notation of section 3. 

Let H(x) = \{x,A{x)), x e G be the reduced Hamiltonian of the left- 
invariant geodesic flow such that A = Al + Ad preserve orthogonal decompo- 
sition G = L + D. Then M = 1(D) = D and Euler-Poincare-Suslov equations 
become: 

£=[£,A D (£)] D , HeD. (22) 

Suppose that L is a Lie algebra of some connected Lie subgroup £. Then 
[L,L] C L, [L,D] C D. Further, suppose that H D {£) = l(€,A D (Q) is an 
invariant of adjoint action of £ on D. Recall that function F : D — > M. is Adz 
invariant if it satisfy equation: 

K,VF(£)]l=0, £eD. (23) 

Let x = £ + rj be orthogonal decomposition of x 6 G, such that £ € D and 
7/ G L. Let A* (a;) = A*(£ + v ) = A D (£). 
Define the function H* : G -> K by: 

#*(*) = ifoA'Oc)) (24) 

Using Ads invariance of Hd(^), we are getting that the Euler-Poincare 
equations for functions H (x) and H*(x) leave the plain D invariant and on D 
coincide with Euler-Poincare-Suslov equations (|22|). For example, the system 

x=[x, VH* («)] = [a;, A* (a;)] (25) 

after projection to the L and D becomes: 

£=[£,A D (£)] + [ri,A D (£)], (26) 

r) = 0. 

Let /i, /i* : T& — > R be the functions obtained by left translations from H 
and _ff * . The nonholonomic geodesic flow has a Hamiltonian restriction property 
with functions h and h*. 

The Hamiltonian function h* has a nice geometrical meaning. Suppose that 
D generate Lie algebra G by commutations. Then the distribution T> is com- 
pletely nonholonomic or bracket generating. The Hamiltonian flow Xh* on T© is 
normal sub-Riemannian geodesic flow of the sub-Riemannian metric induced by 
restriction of the given left-invariant metric to V (for more details see O, M ) . 

We can summarize previous considerations in the following theorem: 
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Theorem 6.1 Suppose that constrained reduced Hamiltonian function Hd(£) = 
2 (£, Arj (£)) is an invariant of adjoint action of the Lie subgroup £ on D and that 
D generate Lie algebra G by commutations. Then on the constrained submani- 
fold T> the following three different problems have the same flow: nonholonomic 
geodesic flow, geodesic flow with Hamiltonian h and sub-Riemannian geodesic 
flow with Hamiltonian h* . 

Recall that equations ( |25|) are completely integrable if they posses a complete 
involutive set of integrals. T = {-Fi, . . . , F^} is complete involutive set of func- 
tions on G if k = i (dim G + rank G) , dFi A ... A dFk 7^ on an open dense 
set IA of G and functions in T are in involution with respect to the Lie-Poisson 
brackets: 

{F h F 3 }(x) = {x, [V-Fj(x), VFi(x)]) = 0, for all x e G. 



If the equations (g5|) are completely integrable then the set W is foliated on in- 
variant (dim G — fc)-dimensional tori with quasi-periodic dynamic (see jl], |2l| ) . 
Also, the system x — Xh* is non-commutatively integrable: T25 is almost ev- 
erywhere foliated on invariant isotropic tori (see p5[). 

Remark 6.1 Integrability of (Eq) do not implies integrability of the Euler- 



Poincare-Suslov equations (22). As an example, it could be happened that 
U n D = 0. If U n D is an open dense set of D, then the Euler-Poincare-Suslov 
equations (E2|) are also integrable. In this case, the reconstruction of the flow 
on the whole phase space V = M. follows from the reconstruction of the flow of 
Eulcr-Poincare equations: the phase space V is foliated on invariant tori with 
quasi-periodic dynamic. 

Example 6.1 Suppose we are given a chain of connected subgroups: 

£ = ©o C ©i C . . . C <S„ = © 
and the corresponding chain of subalgebras: 

L = G c G x c . . . c G n = G. (27) 

Let Gi = Gj_i + Vi be the orthogonal decompositions. Then D = V\ + . . . + V n . 
If we define the left-invariant metric by A of the form: 

A = A L + A D , A D = si-Idv 1 +--- + s n -Id Vn , a,- > 0, i = 1, ...,n, (28) 

then the function Hjj(t;) will be an invariant of adjoint action of £ on D. So 
nonholonomic geodesic flow has a Hamiltonian restriction property. 

Suppose that either (Gi, Gj-i) is a symmetric pair or Vi is a subalgebra of G 
for alH = 1, . . . , n. Let J-q be arbitrary complete commutative set of functions 
on the Lie algebra L (see constructions in uW) lifted to the functions on G. 
Similarly, if Vi is a Lie subalgebra, let T% be a complete commutative set of 
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functions on Vi lifted to the functions on G. Otherwise (i.e., if (Gi,G,_i) is a 
symmetric pair) we take Ti to be given by: 

T t = {f:G^R, f(x)= P (x Gz _ 1 +Xx Vl ), A G R, p6l(G ( )}, 

where X(Gj) is the algebra of Ad@ i invariant polynomials on Gi. 

Then it follows from Mikityuk's results (see JL4J ) that Euler-Poincare equa- 
tions (|25|) will be completely integrable. The complete commutative set of in- 
tegrals is T = Tq + T\ + . . . + T n , 

For example, let us consider the chains: 

L = so(k) C so(k + 1) C . . . C so(n), 

L = u(k) C u(k) + m(1) C . . . C u(n - 1) C u(n - 1) + u(l) C u(n), 

L = sp(k) C sp(k) + sp(l) C . . . C sp(n — 1) C sp(n — 1) + sp(l) C sp(n). 

Corresponding distributions are completely nonholonomic. Whence we get 
integrability of sub-Riemannian geodesic flows on Lie groups SO(n), U(ri), 
Sp(n) with left-invariant sub-Riemannian metrics defined above. 

Example 6.2 Let the Lie subgroup £ be commutative. Then the Lie algebra 
L is contained in some maximal commutative subalgebra K C G. Let D, U 
be the orthogonal complements of L in G and K respectively. Furter, let D 
generate G by commutations. 

Let a be a regular element of K . Then we have that K = {rj G G, [i], a] =0}. 
Let b belongs to K and let R : K -^ K be symmetric operator which preserve 
decomposition K = L + U. By <p a ,b,R denote operator (so called sectional 
operator pl| ) defined with respect to the orthogonal decomposition G = K + 
[a,K\: 

tPa,b,R\K = R, fa,b,B,\[a,K] = ac ^a aa ^b- 

For compact groups, among sectional operators there are positive definite. 
Take such (p a ,b,R- Let H(x) = h{x, f a ,b,R{ x ))- It can be proved that H D {£) = 
!(£> Vo,6,-r|-d(£)) i s an invariant of adjoint action of the Lie subgroup £ on D. 
Thus, nonholonomic geodesic flow has a Hamiltonian restriction property. 

The function H* is of the form: H*(x) = ^(x,(p a ,b,R*(x)), where R* : K — > 
K has a kernel equal to L. For general a G K, the Euler-Poincare equations (E5J) 
are completely integrable both on G and on D. The integrals can be obtained 
by Mishchenko-Fomenko method, based on the shifting of argument of invariant 
polynomials (for details see pl|). Thus corresponding nonholonomic and sub- 
Riemannian geodesic flows are integrable. 
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